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Abstract 

in ■ 

Casimir force encodes the structure of the field modes as vacuum fluctuations and so it is sensitive 
to the extra dimensions of brane worlds. Now, in flat spacetimes of arbitrary dimension the two 

> \ 

standard approaches to the Casimir force, Green's function and zeta function, yield the same result, 

but for brane world models this was only assumed. In this work we show both approaches yield 

| the same Casimir force in the case of Universal Extra Dimensions and Randall-Sundrum scenarios 

\q ■ with one and two branes added by p compact dimensions. Essentially, the details of the mode 

OO ! 

{S| ■ eigenfunctions that enter the Casimir force in the Green's function approach get removed due to 

^: 

\ their orthogonality relations with a measure involving the right hyper- volume of the plates and this 

leaves just the contribution coming from the Zeta function approach. The present analysis corrects 
previous results showing a difference between the two approaches for the single brane Randall- 

^ • 

Sundrum; this was due to an erroneous hyper-volume of the plates introduced by the authors when 
using the Green's function. For all the models we discuss here, the resulting Casimir force can be 
neatly expressed in terms of two four dimensional Casimir force contributions: one for the massless 
mode and the other for a tower of massive modes associated with the extra dimensions. 
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I. INTRODUCTION 



Historically the idea to consider our observable 4D universe as a subspace of a higher 
dimensional spacetime has a long tradition that started with the works of G. Nordstrom 
[lj], T. Kaluza Q and O. Klein |3| (see e.g. ^| and references therein). Nowadays there 
are two broad approaches one typically takes to address the possible consequences of extra 
dimensions in 4D physics. The top-down approach starts either from a fundamental theory 
or a low energy limit of it, for instance M/string-theory or supergravity [5| and upon com- 
pactification of the extra dimensions one hopes to find an effective theory in 4D containing 
as much of the physics we know (see e.g. [6fl and references therein). In this approach one 
favors the properties of the compactification manifold and upon the requirement that the 
compactification be performed in a consistent way, one tracks the physical consequences 
that the geometry of the internal manifold has on the resulting lower dimensional theory, 
including, for instance, the gauge group and the matter content. However, in this approach 
we are unable to select the lower dimensional theory in a unique fashion. The Standard 
Model hopefully would correspond to a particular internal space or vacuum configuration 
chosen by nature by some still unknown mechanism (see e.g. [?] and references therein). 

In contrast the bottom-up approach relies on "model building" , where the requirements 
of having the low energy spectrum and interactions of the known 4D physics put restrictions 
on properties such as the types of singularities, curvature, symmetries, etc., supported by 
the internal space. The constraints are powerful because they hold for a large class of models 
without having to fully specify the compactification details. Of course they are only neces- 
sary conditions, nevertheless they serve as a useful guide in the search for realistic models 
before a complete theory/model can be explicitly constructed. In this approach one looks at 
the different well known physical phenomena and their corresponding experimental confirma- 
tions, and then, by requiring agreement between the contributions of the extra dimensions to 
the 4D physics and the experimental errors, one gets bounds to the higher dimensional free 
parameters. This information forms the core of the necessary knowledge for model building. 



lowing this approach most attention has been devoted to high energy physics (see e.g. 
9) and references therein) and cosmology (see e.g. 10Nl3| and references therein). More 



recently the possibility to obtain information from models with extra dimensions studyin. 



low energy physical phenomena such as the Casimir effect has also been addressed 14j-|34|. 
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The interest in the Casimir force is twofold. Firstly, the force between neutral perfect con- 



ducting plates predicted by H.B.G. Casimir 



35], is experimentally well established |36|-|40|, 



and nowadays the increasing accuracy reached in its determination, makes us think that 
constraining model parameters in this way is at the least complementary to those based on 
high energy experiments. Secondly, its theoretical analysis involves two aspects naturally 
appearing in the study of models with extra dimensions, namely the mode structure of mat- 
ter fields and the submillimeter length scale, of order 1 /im, at which the force becomes 
noticeable and for which some extra dimensional models haven been conjectured to produce 
observable effects. 

In this paper we follow the model building approach to determine the Casimir force for 
a massless scalar field between two parallel plates. This situation mimics the actual exper- 
imental setup where the electromagnetic rather than a massless scalar field is considered. 
We model the plates as codimension one hyper-surfaces in the extra-dimensional space-time, 
therefore what one really obtains is the force per unit of hyper-volume of the plate, as it 
was established long ago for hyper-dimensional Minkowski space-time 41]. In this way the 
extra dimensions yield corrections to the usual 4D Casimir force. Remarkably the resulting 
force can be expressed as the sum of two types of contributions: one that is given by the 
zero mode, thus producing the standard 4D Casimir force for a massless scalar field, and 
the other one that includes the addition of 4D Casimir forces corresponding to the massive 
modes. 

The present work is aimed at showing that both Green's function and Zeta function 
techniques (see e.g. 42] and references therein) yield the same Casimir force for some 
typical extra dimensional scenarios. In particular it corrects a pre vious difference between 
the Casimir force for one-brane Randall-Sundrum models [4i 
method 



16 



44] using the zeta function 



19] and the one obtained using Green's function approach [171.118]. Such difference 



was originated by an erroneous hyper-volume factor for the plates considered in the setting 



in 



17 



18|. 



For the sake of clarity we first study the case of Universal Extra-Dimensions in 5D. This 
corresponds to 4D Minkowski extended by an spatial compact extra dimension attached to 
each of its points. The topology of the extra dimension is an orbifold S l /Z 2 . The Casimir 
effect in this geometry was studied using the zeta function regularization method in [ijj] 
but here we present the corresponding Green's function analysis. The second model we 



3 



shall consider is the so called Randall- Sundrum lip model (RSII-p). These have a single 
(3 +p)-brane 45]. For this model the Casimir effect was computed in [tgI . [19I using the zeta 



function method and in 



17 



1* 



using the Green's function method. Finally we shall consider 



the Randall- Sundrum Ip model (RSI-p). These are defined by two (3 + p)-branes. In this 



case the Casimir force was studied in 



19| using the zeta function technique. To the best of 



our knowledge, an analogous study is missing applying the Green's function approach and 
in this paper we fill in this gap. We shall conclude that for all the above extra dimensional 
models the Casimir force obtained by either of the approaches: zeta function regularization 
or Green's function, is the same. 

The structure of the paper is as follows. Since it will be used frequently in Section [III we 
briefly recall the analysis of the Casimir force in 4D Minkowski space-time for a massive scalar 
field whereas its extension to d + 1 Minkowski space-time is summarized in the Appendix. 
In section II III we discuss the Universal Extra Dimension model. Section [IV] is devoted to 
RSII-p whereas Section [V] deals with RSI-p. Finally, Section [Vl] contains the discussion of 
our results. Unless otherwise stated we use units in which h = c = 1. 



II. SCALAR FIELD IN 4D MINKOWSKI SPACETIME 



To make use of it in the sequel we briefly review the analysis of the Casimir force for a 



massive scalar field in 4D Minkowski spacetime 46]. We start by computing the dispersion 
relation and then determine the Casimir force by the two approaches: Zeta function and 
Green's function. 

Let the scalar field to have mass \i and subject to Dirichlet boundary conditions at the 
planes z = 0,1 . The starting point of the analysis is the Klein-Gordon's equation 

(□ + /i 2 )0 = o, n = d tt -A, (i) 

with A the Laplacian in IB? and // the mass of the scalar field, x represents a spacetime 
coordinate with components (t, x\, x%, z), the last three being spatial and Cartesian. The 
4D Minkowski metric is rj^ = diag{l, —1,-1,-1}. By separating the dependence of the 
field in Cartesian coordinates x = {x\, x 2 , z) as <p(xi, z, t) = Xi{xi)T(z) e tult , i = 1, 2, we have 
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the equivalent set of eigenvalue equations 

-duXi(xi) = fsfxifa), i = 1, 2 , (2) 
-d zz T(z) = k 2 3 T(z), hi := cu 2 - (/i 2 + k\ + k\) . (3) 

Here k 2 , i — 1, 2, and are separation constants or eigenvalues. The physical plates are 2D 
surfaces described by coordinates x\ and xi so that Xi in © can be subject to free boundary 
conditions. Coordinate z is transverse to the plates so T in ([3]) will be subject to Dirichlet 
boundary conditions at z — 0, Z. The corresponding eigenfunctions are 

Xi(xi) = -^=e ikiXi , h eM,i = 1,2, (4) 
v 2ir 



T N (z) = A/ysin^p, A = l,2,..., (5) 



and the resulting dispersion relation is 



co 2 = /i 2 + k\ + A; 2 + (6) 



A. Zeta function approach 

To compute the Casimir force one can compute the Casimir energy between the plates 
Epiatesi by summing up the zero-point energy per unit area hu/2. There are two ingredients 
required to follow this strategy: the dispersion relations and the modes structure (which can 
be continuous, discrete or an admixture). It turns out that E p i ates contains a linear term 
in the separation I between planes which gives rise to a constant Casimir force. This term 
can be canceled by addition of a constant to the Hamiltonian density or by considering 
the energy E in the absence of the plates which means k% G M and T^(z) = -j=e tk3Z . 



12-k 

Adopting the second option, the resulting finite expression for the Casimir energy per unit 
area of the plate is 



c i \ Epl a t e s — Eq 1 T-r f°° dki I ST^ , \ , f°° dk 3 

em = — T2 — = 2 n i l - 1 y ^ ^ 

i=l,2 J ~°° \N=1 J ~°° 



k 2 ,k 3 (^) » ( 7 ) 



where 



A 2 7T 2 

u kxMjN {n) = \l fi 2 + kj + k 2 + —p-, (8) 



00 = y/f>? + k\ + k\ + k 2 , (9) 



and L 2 is the area of a square shaped piece of the plates at z = 0, 1. Notice the extra factor 
of I in the second term of (jTJ); it comes from the fact that Eq is the energy in the whole 
volume delimitated by z = 0, 1, in the transverse direction, whereas £ is the energy per unit 
area L 2 . We are denoting explicitly the dependence of the density energy E±d on the mass 
/i, to stress the fact that the 4D scalar field is massive. 

We perform explicitly the integrals in the appendix IVHI obtaining 

2 1 00 1 

^) = -^-- l Y,jp K -^ m ^ ( 10 ) 

JV=1 

where K is the modified Bessel function of second type. The Casimir force is obtained 
from the Casimir energy simply deriving with respect to the separation between the plates: 
F 4 d = —dSio/dl, thus 

(11) 

where a property of the derivative of the Bessel function has been used. In general this 
expression can not be simplified further and usually people computes it numerically for a 
given value of the mass \i. For the massless = 0, however, the expression can be 

simplified to yield 

M0) = ~, (12) 
in which the appropriate approximation for small argument of the Bessel functions has been 
used and then the identification of a zeta function allows to evaluate the result. 



fAD{lA 



1 1 9 00 1 

. N=l 



I ^ N 

N=l 



B. Green's function approach 

In the Green's function approach, once we are armed with the eigenfunctions (jlj) and fl5]), 
we can express the Green's function G^d for the problem (□ + / u 2 )G , 4 £) (x, x') = —S(x — x'), 
subject to Dirichlet boundary conditions at z — 0, 1, as 

G iD (x,x') = II / dk iX *(xi)Xi(xd [ pe-^-^g&z'), (13) 

i=l,2"' 

a (^') = ± T ^™f , (14) 

JV=1 I 2 ^3 

with * denoting complex conjugation and g the so called reduced Green's function. Notice 
that this expression of the Green function is valid only in the region between the plates. 
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Now given the relation between the vacuum expectation value of the time ordered product 
of fields and the Green's function, (T[(j)(x)(j)(x')]) = ^G^d^x, x'), the force per unit area on 
either plate can be obtained from the vacuum expectation value of the energy-momentum 
tensor T up = d v (j)d p (j) — r\ vp L with C = \d u (pd u (f) — \p?^) 2 . Since upon integration over all 
space the term r] up C does not contribute by virtue of the Klein-Gordon's equation ([T]) one 
gets 

fM = (T™)= ]imd z d z> G AD (x,x') (15) 

x'-^x 2= o/ 



dkiX*( x i)Xi(xi) / — lim d z d z ,g(z,z') (16) 

J ATT z'^z z=0 j 



In the coincident limit Xi( x i)Xi( x i) — ^ an d so the dependence on x^i = 1,2, drops out; 
which should have been expected from the translational invariance of the parallel plates 
configuration along the Xi,i = 1,2, directions. 

Combining ([5]) together with ( Fl4l) allows us to obtain the explicit form of g(z, z'), which 
upon substitution in ( Fl6|) and after the change of variables: u — > i(, and k\ + k\ + ^ 2 — > p 2 
produces 

1 f dk x dk 2 f d£ A , , 

This allows us to read (fc«,0 as three-dimensional Cartesian coordinates. This is not the 
final answer because so far we have only considered the force to the left of z = I or to 
the right of z = 0; actually the integral ( |T71) diverges. We also have to include the flux of 
momentum for instance to the right of z — I. We shall not elaborate on this issue, for our 
purpose it is enough to mention that the normal-normal component of the stress tensor at 



z — I is: (T™*) = iyj p 2 + /i 2 /2 46_|. The net forces at z = I produces finally 

1 fdkidh fd( (2^k rj- — A , 1a . 

f«M = 2j w/^i'g 1 ^ + ^ + ^ + V7T ^J- (18) 

We compute this integral explicitly in the appendix IVHl obtaining again the expression ( TTTT) 
for the 4D Casimir force. We could use the easier argument based on the discontinuity of 
the derivative of the reduced Green's function g above. Of course the same results from 
the discontinuity of the zz component of the energy momentum tensor on either plate at 
z = 0,l. 
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Notice that in order to obtain the force, the Green's function method includes a "bit" 
more of information than the previously used in the zeta function regularization method. 
We have used again the eigenvalues through the dispersion relations but we have also used 
explicitly the eigenf unctions and not only the mode structure. However upon integration 
of the modes Xii * = 1, 2, the real input is again only the modes structure as in the zeta 
regularization method loosing the "bit" of extra information. This property of the Green's 
method is the one that makes it equivalent to the zeta function regularization method. 
For instance in the models with extra dimensions, as we will discuss, the eigenf unctions 
depending on the extra coordinates give us information about the localization of the field 
modes, however, since at the end these eigenfunctions are integrated out, the information on 
the localization is in some sense "lost" . This is in agreement with our concept of consistent 
compactifications, for which the extra dimensional coordinates must disappear explicitly 
(see for instance 47| and references therein). 



It is straightforward to generalize the result of the 4D Casimir force to a (d + 2)D 
Minkowski spacetime (see appendix IVIII) . In this case the scalar field is bound by hyper- 
planes of d dimensions and the force per unit <i-dimensional volume between the hyperplanes 



is given by 



4l| 

" 3 ^ 1 T ^ /oAr; . 1 



f(M)o(ji) = -2 (£ 



d+2 



(19) 



For the forthcoming analysis it is convenient to notice that this expression has the following 
limit values 

lim WM = -j^ElT p±*) C(d + 2), (20) 

and 

lim f {d+2)D {y) -> 0. (21) 

Once we have reviewed the way in which each method gives origin to the Casimir force, let 
us continue with the extra dimensions models. 
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III. UNIVERSAL EXTRA DIMENSIONS 



The model 



The Casimir force in a Universal Extra Dimension (UXD) scenario 48| was considered 



in 



14j for the case of a massless scalar field to probe the possible existence and size of 
an additional spatial dimension which is compactified on a S 1 /Z 2 orbifold. This geometry 
restricts the possible vacuum fluctuations of the scalar field to have a wave vector along the 
extra dimension of the form k n = n/R, with k n being the wave vectors in the direction of 
the universal extra dimension and R the radius of S 1 . 

Let us start with the 5D action for a massive scalar field 

S=^Jd 4 xj^ Rd9^(g a Pd a <$>(x,9)dp$(x,9) - m 2 5 <$> 2 (x,9)) . (22) 

Here x a = (x^, R9) are the coordinates with a = (/i, 4) and \x = 0, . . . , 3 are the indexes of 
our 4D spacetime. m 5 is the mass of the 5D field. ttR is the size of the extra dimensions 
and g afi is the inverse of the metric defined by the interval 

dsl = r]^dx^dx u - R 2 d9 2 . (23) 
In this metric, the 5D Klein-Gordon equation reads 

□ 4 $--^<9 e 2 $ + m^ = 0, (24) 

which separates through <&(x,R9) = 4>{x)ip{9) into 

{d 2 + m 2 R 2 )m = 0, (25) 
□ 4 + (m 2 5 + m 2 ) <j)(x) = 0. (26) 

When the extra dimension is Kaluza-Klein type, the only condition on the fields is: $(x, 9) = 
<&(x, 9 + 27r), which allows for a Fourier expansion taking the form 

i 00 nr 

$(x,0) = ^ (0) (i)+yj — \<f) {n \x)cos{n9) + x (n) (x)sm(n9)} . (27) 
VttR ^ V ttR 

When the extra dimension is instead an orbifold S 1 /Z2, there is an additional parity condi- 
tion on the 5D scalar field: <&(x,9) = ±$(x, — 9). In this case it is clear that the modes of 
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the scalar field have definite parity, they are either even, denoted by (+), or odd, denoted 
by (-). Explicitly, they are 

4>(x,6) + = -= 4>W(x) +^ J— 4^{x) cos(n0), (28) 

v n=l 

<j>{x,9). = J2\I^X {n) (x)sm(n9). (29) 

n=l 

In models of extra dimensions, the zero mode is associated with the lower dimensional 
physics, and, as a consequence, we need to consider the even modes if we want to repro- 
duce the lower dimensional physics. Therefore we impose the additional parity condition 
$(x, —6) = +<&(x,9) on the 5D scalar field which leave us with the set of 4D scalars field 
{(j)(°\x), (fi( n \x)}. These fields satisfy the effective 4D equations 

D 4 (j) {0 \x) +m^(°)(x) = 0, (30) 
U 4 <p {n \x) + [m\ + ^ 4>M(x) = 0. (31) 

We interpret to the zero mode <fy^>{x) as a 4D massive scalar field of mass and to each 
mode of the Kaluza-Klein tower 0^ n ^(x) as 4D massive scalar fields of mass 777,4, given by 



n 2 



m 4 =\/ml + —, nGNU{0}. (32) 



The modes for the coordinate 9 are simply 



^(6) = -^= and Vn(0) = \/-^cos(n0). (33) 

V7T-R V ttH 

This means that the zero mode field is constant along the extra dimension whereas the 
massive modes are distributed harmonically in that direction. 



B. Zeta function approach 

In [14J] the zeta function method was used to compute the Casimir force associated to a 5D 
massless scalar field. Here we consider a 5D massive scalar field. As we have mentioned, in 
this formalism the relevant quantities are the frequency of the vacuum fluctuations and the 
modes structure, so if we decompose the 4D fields (i.e. the zero mode and the Kaluza-Klein 
tower) in the same way as in equations (|2"1) to (131), the energy per unit 3- volume (L 2 x irR) 
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of the hyperplanes is 



£uxD(m 5 ) = i Yl I T 1 f £ Uk u N,n(m4)-l [ ^u^fe^N , (34) 



where 



Wfc s ,Af,n (m 4 ) = \ kf + k?, + 



7TN 



n- 



n- 



Uki,k 3 ,n (m 4 ) = \/ fc? + fcf + fcjj + + m 5- 



(35) 
(36) 



and the mass 777.4 is given by fl32l) . 

In order to get £uxd we do not have to compute anything, we already have the answer. 
Notice that we can rewrite Eq. (|34p as follows 

00 

(37) 



£uxD{m 5 ) = £ 4 £>(m 4 | n=0 ) + ^£iD (m 4 ) 



n=l 



where £40 (^5) is given by (jTJ) and whose analytical expression after integration is ffTO]) . 
Deriving this expression with respect to the separation between hyperplanes we finally get 



fuxD(m 5 ) = f 4D (m 5 ) + ^ f±D (m 4 ) 



(38) 



n=l 



which reads 
fuxD(m 5 ) 



mi. 



■in' 



3 



00 , „ 00 

£ —K&Nlms) + E ^1(2^775) 

JV=1 



7Y=1 



n=l v 7 L iV=l \ / 



+ 

For the massless case (7775 = 0) [14], one obtains 

2 1 

7T 1 



2^/777! + ^ 



AT=l 



77/ 

R 2 



fuxD(m 5 = 0) 



480 / 4 
480Z 1 



77 



+ £/">(t;) (39) 

n=l 

21 f 4 4 r A ^ (^-) +-t 1 k 1 (2n£) 

8vr 2 ^ R 2 I 2 ^ N 2 2 V RJ IR^ N 1 V i?/ 

n=l L AT=1 AT=1 

In the setting with only one extra dimension, a good agreement with the data can only 
be obtained if the radius of such dimension is smaller than R <10nm [14]. This bound is 
weaker than others obtained from high energy physics which are around 10 _9 nm. 
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C. Green's function approach 



Consider now the 5D Green's function expressed in terms of the eigenf unctions (|33|) 

G 5D (x, 6, x', 6') = J2 M0)M#)G4D [x, x'- ml + ^y (40) 

where G±d is the 4D Green's function given by (|T3|) . In terms of it, the Casimir force between 
the hyperplanes is 

fuxo = ^jd^ [ Rde^9 [{TZ)l=i - (T™%=i\ , (41) 



where 



{T in/ou %=i = La^G^ix, 9, x', 6') \ Xi ^, e ^. (42) 



Just as with the zeta function method, we do not have to compute that much to get the 
answer. Notice that the 5D Green's function can be rewritten in terms of the 4D one in the 

way 

1pn{VWn{P ')^4D {X,X;m~ 5 + — 
n=l " 

and the expectation value (H2|) can be rewritten in terms of the expectation values in 4D 



G 5D (x, 9, x', 6') = Md)M0')G4D {x, x'- ml) + M0)M0')G4D (x, x'; ml + ^y (43) 



00 

n=l 



Substituting this expression in (j4ip we can rewrite the force in terms of the Casimir force 
for scalars fields in 4D 

fux D (m 5 ) = / 4D (m 5 ) jf RdB$(6) + (m 2 5 + ^ jf Rdd^). (45) 

Because both integrals in 9 are equal to 1, we get exactly the expression f )38|) . We conclude 
that in the case of one Universal Extra Dimension, the effective Casimir forces obtained by 
both methods coincide. Notice that he force we have computed is the force per unit volume, 
i.e. the force per unit area (L 2 ) of the plates and per unit length in the extra dimension. 
Physically what we have is a couple of 3D plates, with one dimension stretching along the 
extra dimension, but both embedded in four spatial dimensions. Such a setting is referred to 
as having plates of codimension one. An extension of this idea to Randall-Sundrum models 
actually holds and is what we show next. 
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IV. RANDALL-SUNDRUM II-p MODELS 



The interest in the Randall-Sundrum II-p models comes from its property of localizing 



not only sea' 
dimensions 



ar and gravity fields but also gauge fields whenever there are p extra compact 



45 



50(. In the case of p = the model only localizes scalar and gravity 



fields. The model corresponds to a (3 + p)-brane with p compact dimensions and positive 
tension k, embedded in a (5 + p) spacetime whose metrics are two patches of anti-de Sitter 
(AdSs +p ) of curvature radius n~ l 

p 



Ub 5+p 



-2 K \y\ 



(46) 



The Casimir force for a massless scalar field in the RSII setup (p = 0) was computed 
in [lfj using the zeta function regularization method, whereas in the Casimir force 
was computed for both a massive and a massless scalar field in the RSII-1 model and then 
generalized to the RSII-p model by means of the Green's function approach in [3] and using 



the zeta function method in 



191 ] . However their results turned out different and seemingly 



depended on the method adopted - a situation clearly unacceptable. Here we will show that 
such difference was originated by an erroneous hyper- volume factor for the plates considered 



in the setting in 



17 



18l | . We shall restrict ourselves to the case of a 



ligher dimensional 



181 ] performing the 



massless scalar field. The interested reader in the massive case can see 
corresponding modifications. A related calculation of the Casimir effect in de Sitter and 
anti-de Sitter braneworlds can be found in 5l| . 



A. The mode structure 

Let us consider the (5 + p)D action for a massless scalar field $ in the RSII-p metric (|4"6"|) 

1 r P 

S=- d 4 x\\ RjdOj dy^/\g~\g aP d a § dp®. (47) 
J i=i 

Here X a = (x^, Rfii, y) where x^ are the coordinates of our 4D spacetime, the p coordinates 
9i are associated to the p compact S' 1, s and y is the noncompact coordinate transverse to 
the brane which is placed at y = 0. The field equation for the scalar field is given by 
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which separates through $(X) = ip(x) Yl^=i ^ > j(^j)i J (y) m t° 

(dl + mlRf) Qtf-) = 0, j = 1, . . . ,p, (49) 
(d*-(4 + p)Ksgn(y)d y + m 2 e 2KM ) ip(y) = 0, (50) 

(n 4 + ml)tp(x) = 0. (51) 

The (p + 1) separation constants with units of mass, mg. and m, correspond to the spectra 
of the modes for the compact and non compact dimensions, respectively. They give rise in 
turn to the effective mass, m 4 , of the 4D modes in ( l5ip through: m\ = Y^j=i m e + m ' 2 '■ 

To find mode solutions to the above equations we shall incorporate three types of bound- 
ary conditions: (a) To implement the presence of the plates in (4 + p)-space we simply set 
(p(z = 0, 1) = 0. The eigenfunctions and eigenvalues for this Dirichlet boundary conditions 
were already discussed in section [Til (b) To match the modes across the brane along the non 
compact dimension we impose ij){y = + ) = ip(y = CT) and d y ^{y = + ) = d y ip{y = CT). (c) 
To account for the compactness of the p dimensions we set Q n .(9j) = ® n A6j + 27r). Hereby 
we obtain explicitly the plates represented by two parallel planes in 3-space but stretching 
along the extra dimensions. 

The allowed modes for the non compact dimension are now a massless zero mode localized 
on the brane 

*■ - (52) 

which satisfies the normalization condition, 

/oo 
dye -(p+*)K\y\^ = l. (53) 
-oo 

The localization comes from the fact that the 4D effective profile of the modes is given by 
ipo = e~( p+2 ) re ' J/ '' 2 ipo which clearly is localized on the brane. The massive modes have the 
form 

m > 0. (54) 



a m J 7 ( — — ) + TO iV 7 



Here J 7 and N 1 are the Bessel and Neumann functions respectively. 7 = and the 
coefficients a m and b m are given by 

A 1 

(55) 
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where 



A, 



^7-1 (*) 



(56) 



Notice that in this case the localization of the massive modes on the brane is better for 
increasing p, since the modes are modulated exponentially in the form e~ pK ^ y ^ 2 . The nor- 
malization condition for the massive modes is 



The modes in 9j are: 




1 



where 



rij = m e .Rj G Z. 



(57) 



(58) 



Therefore the contributions of the extra compact dimensions to 7714, are given in terms of 

(59) 



ml = nVm 



v 2 
3=1 



B. Zeta function approach 

In analogy with the cases studied above and due to the fact that the modes behave 
differently for the zero mode (1521) and for the KK modes (I53|) . the energy density per unit 
of (p + 3) volume (l 2 x n?=i( 2 7r#j) x (p+3) K J ^ or the sca l ar held is 



£ 



RSIIp 



^n/ltE e 

+ 2 11 / 2W 72. L (^) - « / 



(60) 



(m 4 ) 



where m 4 is given by fl59|) . {n} denotes the set {n 1; n 2 , . . . , n p \rii G Z, . . . , n p G Z} and the 
dispersion relations are 



^ki,N,ni,m ( m 4) = 



fc^ -|- ~t~ 



TriV 



3=1 n 3 



(61) 



and 



(m 4 ) 



P ™ 2 

kf + k% + k 2 + + 



m 2 . 



(62) 
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It is important to stress that in this case each one of the different p sums in rij goes 
from — oo to oo and not as in the UXD case from to oo. The reason is that here we are 
considering the Kaluza-Klein tower associated to S 1 whereas in the UXD case the tower is 
due to the orbifold S 1 /Z 2 . In fact the Kaluza-Klein tower of S l is two copies the Kaluza- 
Klein tower of S 1 /Z 2 . 

Regarding the integration on the continuous massive modes to it is possible to take 
advantage of the following trick 

f°° 1 f°° f°° dm 

/ dmf(m) = - dmf{m) = 71 —/(to), (63) 

J0 2 J -oo J -oo to 

which is valid whenever the function / be even: /(—to) = /(to). Due to the fact that the 
frequency a; (to) satisfies this condition, one can consider the integration on to at the same 
footing that the integrals on ki, i — 1, 2. As a consequence after integration, the result is 
simply 

WO) = g ( + > ( ) , (64, 

where £5.0 (/x) can be obtained from ( 11031) and whose analytical expression after integration 
can be obtained from (11051) . Deriving this expression with respect to the separation between 
hyperplanes one gets 



fnsiip(0) — ^2 ( ^ 4D 



Interpretation of this result is straightforward, the first term corresponds to an infinite sum 
of 4D Casimir forces, one of the terms corresponds to the Casimir force due to a massless 
scalar field (the one corresponding to the zero mode of all the p compact extra dimensions), 
plus an infinite sum of 4D Casimir forces corresponding to massive scalar fields (where the 
mass corresponds to all different combinations where there is at least a non zero mode), 
the second term comes from the non zero modes of the non compact extra dimensions and 
corresponds to a sum of 5D Casimir forces. Some examples previously discussed in the 
literature are: 
• Case p = 0. 

In this case we simply have [if]] 

W(0) = /id(0) + -/ 5 d(0). (66) 
16 



v 1 
n 



At this point you can think units do not match, but they do because is a force per unit 
area and /s is a force per unit volume. Explicitly we have 

fRSIl(0) - 



TV 2 1 



45 1 
480 Z 4 V 4^3 C( Q 



(67) 



Case p = 1. 
In this case we have 



12| 



n 
>7? 



n=0 



(68) 



Wi(0) = J2 (i) + ^ ( 

n=— oo 

oo 

= Ad(0) + IfMo) + 2^ (Ad (£) + -J,d ( 
C. Green's function approach 



Lets now apply the Green's function method to the RSIIp models. As we have discussed 
the force between the plates is obtained by integrating over coordinates "lateral" to the 
plates. In this case: x±,y,9j due to the fact that the normal- normal component of vacuum 
energy momentum tensor in 3+1+p spatial dimensions has physical units of force per unit 
of "volume" of 2+1+p space: 



/•A poo i 

/ dx± / dyJ\g p i ate \ 


P r2w 

n / 










h J-oo v 


b=i j ° 






z=l 


z=l- 



(69) 

where A is the area of the planes forming the plates in 3-space and g p i a tes is the induced 
metric on the physical plate located a z = I. Since the physical plate is a surface of (p + 2)D, 
the jVflbk 



in 
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contributes an exponential of — n\y\ (p + 2). Instead of this factor the exponential 
18] contained an erroneous power given by —K\y\{p + 3). 
The vacuum expectation values and the Green's function are related to the normal-normal 
components of the vacuum energy momentum tensor through 



(T in/ou t) 



dzdz'G^^x, y, 6; x', y', 9' 



z=l 



2i 



(70) 



x i —*x'. , z-^fz'=l. 0,-— >(?'. 
X' ' ■> J 



As in the previous cases we can rewrite the (5 + p)D Green's function in terms of the 4D 
Green's function 

v 

(x,y,6;x',y',6') = Y,U%^ Q -M)My)My')GM 

X, X j ?7i4| m= o) 

{"} i =1 

+ E / ^n©; j (^)0n J (^)^ m (y)^ m (l/ / )G'4D(x,x / ;m 4 ),(71) 
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where m 4 is given by (|59|) . Introducing this expression in (170|) and then in (|69|) we obtain 
/ = £ / 4 D(m 4 | m=0 ) r^e- 2 ^^ +2 Vo 2 + E f—fMm,) f°° dye'^^Uy). 

M • / -° {n} 7 ^ • / -°° 

(72) 

But by virtue of the relations f lo^j) and ( 15"7"|) . the dependence on the y coordinate drops out 
completely obtaining 

f = Y1 /4D(m 4 | m =o) + [ — /4c(m 4 ). (73) 

{n} {n} J ° K 

This expression coincides exactly with the one obtained by the zeta function method f )65|) 
once one performs the trick (|63|) . 

Before ending this section, some comments are in order. Notice that there is no factor 
depending on the number of compact dimensions p in front of the effective 4D Casimir 



force, as it was reported in 



17 



18| . Also, because the integration in the non-compact extra 



dimension has been carried out, it is not longer necessary to evaluate the eigenfunctions 
^miy) on the brane. In fact the contribution of the whole size of the plates in the y direction 
has been considered. 



V. RANDALL-SUNDRUM Ip MODELS 



A. Mode structure 



As a final example we discuss the case of a bulk scalar field in a RSIp model. In this case 
the metric is given again by (Ho]) , but this time the coordinate y is compact (0 < y < irr). 
Thus the setup allows two (3+p)-branes to lie, respectively, at y = 0, 7rr. The Casimir effect 
for a massless scalar field in this model was computed using the zeta function method for 
p = in [l6| and for arbitrary p in |19l| . We address this problem here for a massive scalar 
field of mass /i. As is well known [52], when the higher dimensional scalar field is massive 
there does not exist a zero mode solution of the equations of motion with simple Neumann or 
Dirichlet boundary conditions. In order to overcome this problem it is necessary to modify 
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the boundary action and include boundary mass terms 



S — S$ + S) 



i OB ranc 



(74) 
(75) 



= \ I d i xdyf[R j de j ^(g MN d M ^d N ^-^ 2 ) , 

J 3=1 

r p 

5 Bra nc = - / d A x dy J] Rj ddj yf=g 2bK [5{y) - 5{y - tit)] $ 2 , (76) 

J 3=1 

where h is a dimensionless constant parametrisation the boundary mass in units of k. The 
S-Br&ne allows to implement various boundary conditions corresponding to different mode's 
localizations in the y direction as we explain now. The resulting field equations are, by 
letting $(x, y, 9j) = <p(x)i/>(y) YT j= i 



• ml nj) = (77) 

(d 2 y - (4 + p)Ksgn(y)d y + m 2 e 2K ^ - /i 2 - bK(5(y) - 8(y - irr))) ip(y) = (78) 

(□ 4 + ml) <f)(x) = (79) 

with mg j , m separations constants so that the effective mass of the scalar field can be read 
as m 2 := mj. + m 2 . Eigenfunctions and eigenvalues of the p coordinates 6j are given by 



( |58i) and we do not elaborate further on them. In the y direction the eigenfunctions are 
accounted for by subjecting them to the modified Neumann boundary conditions 



^ - bKsgn(y)ifj 



0. 



(80) 



y=0,7rr 



If we write the higher dimensional mass in units of k, i.e. fi 2 = cik 2 , the above equations de- 
pend on the two arbitrary mass parameters: a and b. For generic values of these parameters 
there are not solutions to the boundary conditions, however if b = a ± 7 where a = -^p and 
7 = Va 2 + a, a zero mode solution exists. Assuming 7 to be real, the only free parameter 
has a range — 00 < b < 00 and using it, the scalar zero mode can be localized anywhere in 
the bulk 



(6-[a-l])w (b-\a-l])KV 
f e 2(b-la-l]) K *r_i\ c ' 



12-KT ' 



|Mtt-l]|* 



1-e 



-216 — \ot — 1] I nirr \ 



,(b-[a-l])Ky 



b — [a — 1] > localized in IR brane 

b — [a — 1] = no localization (81) 
b — [a — 1] < localized in UV brane 
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which satisfies the normalization condition 



I* TIT 

/ dye-^\y\^{y) = l. (82) 
Jo 



These results generalize the case p = (=>- a — 2) [521 ]. 
As for the massive modes we have 



Mv) = e ° KV [ c i J 7 (7^) + c 2 F 7 (™e**)] . (83) 



where are arbitrary constants and J 7 , F 7 are Bessel's functions of order 7. Imposing 
the boundary conditions in the low energy regime m « K with large r: kt » 1 allows to 
obtain the approximated Kaluza-Klein mass spectrum 

/ "" f n+ Z_^\ 7re -' e " T , n = l,2,... . (84) 



k V 2 4^ 
Thus we have that the 4D mass is given by 



TO = < 2 J o (85) 



The constants in ( 1831 are chosen in such a way that the orthogonality relations are 

dye^ 2+ ^Mv)M(y) = <W (86) 



B. Zeta function approach 



Now the energy density per unit of (p + 3) volume of the plate, which is given by 
(l 2 x ri;=i(27r#j) x ^^(1 - e -^ +3 ) K7rr )) > takes the form 



Si*/, - i n /^frE (f>^ (e^ 

i^* 7 " 00 { nj } \JV=1 \i=l J 



/ i=l,2' / -'* ^ n=l { „ j} \ N =1 \i=l ^ / 

/°° dk ( p n 2 \\ 

^^uks^rnl^^+mljj , (87) 
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where {rij} denotes the set {ni, n 2 , . . . , n p \ni 6 Z, . . . , n p 6 Z}, 



p n 2 



\ 



k^ -\- &2 "I - 



vriV 



(88) 



and 



P 2 



Next, upon integrating, we have 



\ 



v i 

Ti 

k\ + e 2 + ki + J2^2+ <■ 

j=i j 



(89) 



P 2 



n=l 



V 2 

^ R 2 



\U 



mi 



(90) 



where S^d is given by ( flOl . Deriving this expression with respect to the separation between 
plates one gets 



Insipid) — E ( f AD 

{Uj} 



\1 



1 / n =l 



P 2 
n^ 



\u 



(91) 



As a particular case we have fi = p = 

/WO) = Ad (0) + f; Ad ( (n + ^ Tre"^ 



(92) 



7T 2 e" 2 ^ <x 



480Z 4 8 



n=l 



3 A i 



J2 W2 K 2 (2Nlln + 



I 2 ^ N 2 

N=l 



7TC 



2 / r 



T E^( M (» + i) 

N=l \ \ / 



7TC 



C. Green's function approach 



Once again the calculation using the Green's function method relies on the fact that the 
(5 + p)D Greens function can be written in terms of the 4D Green's function via 

p 

G { r 0+p)D (x, y, 6; x', y', 9') = EII 6> n j (^') e> n j (9' j )il) (y)ipo(y')G4 D (x, ar'; m 4 | m =o) 

oo p 

+ Y,Y,U%^ e ^M)My)My')G 4D (x,x'-,m,), (93) 

{n^} n=l j=l 
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where is given by (|85|) . Introducing this expression in (170|) and then in (|69|) we obtain 

00 /-7rr 



/=£/ 

{n 3 } 



4Dl^4|m=oJ / dye 





{n 3 } n=l ^ 



(94) 



And now, by virtue of the orthogonality relations, Eqs. ( 1821) and (1861) . the dependence on 
the y coordinate drops out to yield 



/ = E /^(^({nj}, ra = 0)) + E E Ao( m 4({ n i}» 

{n} {n} n=l 



(95) 



This expression coincides exactly with the one obtained by the zeta function regularization 
method, Eq. §2). 

VI. DISCUSSION 

The old idea that our world is embedded in a spacetime with dimension higher than four 
has reemerged in brane world models which have revealed windows to look for deviations 
from standard physics mostly in high energy physics [§, 0] and cosmology (e.g. 
Nevertheless low energy tests may also provide some insi| 
dimensions including in particular the Casimir force fl4 - 



10 



13j). 



it into possible imprints of extra 

n 

-I34j. Such force is sensitive to 



the mode structure of the field which in turn depends on the features of the background 
spacetime and bounds can be set for the values of the parameters of given brane world 
models which produce Casimir forces deviating from known data beyond the corresponding 
uncertainties. 

To determine the Casimir force one can made use of either of two well known approaches: 
Green's function and Zeta function. In the case of flat spacetimes both approaches yield the 
same result (See eg. 46[), however, for brane worlds one usually assumes this is the case. In 
this work we have actually shown Green's function and Zeta function yield the same Casimir 
force for the case of Universal Extra Dimensions and Randall-Sundrum models with one or 
two branes added by p compact dimensions. These results correct in particular an erroneous 
difference between the Casimir force obtained by Green's function technique 13, ll8j and the 



one obtained from zeta function 



m 



191 ] for a massless scalar field in the case of a single brane 



Randall-Sundrum scenario added by p compact dimensions. The origin of the difference in 
this case was due to an incorrect hyper-volume for the plates subject to the Casimir force 
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in 113, |1S | . 

The coincidence of the above two approaches to the Casimir force can be understood 
as follows. Although the Green's function technique involves further details of the mode 
decomposition of the corresponding fields, it is due to their orthogonality relations which 
involve the correct hyper-volume factor of the plates that one can literally eliminate the 
mode eigenfunctions from the Casimir force. This is neatly seen in the case of UXD, Eq. 
(j45p . This holds similarly for the cases of RSII-p and RSI-p, as can be seen from Eqs. f J72|) 
and (I94"|) . respectively. 

Hence given the equivalence of the zeta function and Green's function to determine the 
Casimir force in Randall- Sundrum and Universal Extra Dimensions models one can conclude 
that localization of the field modes does not play a role as far as the Casimir force is 
concerned. This is so due to the fact that zeta function is actually insensitive to the form of 
the mode eigenfunctions but only to dispersion relations. On the side of the Green's function 
approach, while built explicitly on an eigenfunction expansion, it looses them by virtue of 
the orthogonality relations they fulfill which just absorbs the hyper- volume of the plates. 
This is reassuring since it has been noticed that localizing all the fields of the Standard 
Model to the brane located at y = irR leads to problems with the phenomenology of proton 
decay, Flavor Changing Neutral Currents (FCNC) effects and neutrino masses 52| . Indeed, 
although originally the p extra compact dimensions were added to the Randall- Sundrum 
models to produce localization of gauge fields, such a feature seems not to be required 
anymore. 

It should be stressed an important pattern has put forward in this work for the resulting 
Casimir force in brane worlds, namely, it can be expressed as the sum of two terms each one 
having the specific form of the 4D Casimir force: the first one containing in particular the 
zero mode defined by the extra dimensions and the second one including the full Kaluza-Klein 



tower of massive modes. This feature simplified importantly the analysis and in particu 



allowed to adopt the full machinery of previous results in d + 2 Minkowski spacetime 41] 



ar 



It would be rather interesting to study other brane world models to test the equivalence 
we have here proved between the zeta function and Green's function approaches in the 



calculation of the Casimir force. Indeed, it shou 
at least for sufficiently symmetric spacetimes [53]. 



d be possible to have a general proof of it 
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VII. APPENDIX: CASIMIR FORCE IN D+2 MINKOWSKI SPACETIME 



In this appendix we compute the Casimir force for a massive scalar field of mass /x in 
(d + 2)-dimensional Minkowski spacetime. We discuss the calculation using the two methods 
we are interested in: the zeta function method and the Green's function method. Although 
for a given mass \x 7^ the force can be computed at the very end only numerically, it is 
possible to give an analytical expression of it for a generic mass. The aim to show this 
computation is twofold: i) to avoid unnecessary repetitions of the calculation throughout 
the paper and ii) to allow us comparison with recent results in the literature. 



A. Zeta function approach 

This computation was originally discussed in 
Our starting point is equation (j7j) 



41] and we adapt it here to our notation. 



n _ Elates ~ Eq 1 -i-i- f°° dki I ST^ , v , f°° dk^ , . 

8 " J} 2 J- X / ^ L u *to#W ~ l 2^*1*,* UO 

„■ — i o v — OO \ at 1 o — OO 



(96) 



i=l,2 " \JV=i 

In terms of the integral 



J ^ 2) "^/(2^^ + a2 ' (97) 



we can rewrite equation (196)) as 



00 / at-2 2 \ 

E = E J ^ -F2- + ^ - 1 W A (98) 

N=l ^ ' 

In order to be general, we shall evaluate the Casimir energy for a scalar field in D = d + 2 
dimensional Minkowski spacetime, between hyperplanes of dimension d. This problem is 



known as the Casimir effect of codimension one 



4l| 



s d+ 2 = Eplates L ; Eo = f>(^+/H (") 



N=l 

Using the Euler representation for the gamma function 



POO 

F(z) =g z e~ 9 H z - l dt, 
Jo 



(100) 



the integral (|97|) can be rewritten employing the Schwinger proper time representation for 
the square root as 

i r d d k r° dt 



dK } 2r(-l/2)7 (27r) d 7 t 1 J 
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Performing the Gaussian integral first and using (11001) again we have 



h(a 2 ) 



I A \ Cl + 1 



t 



d + l 



a 



d+l 



(102) 



where we have used the value r (— |) = —2y/ir. Substituting this result in (|99|) we obtain 



S, 



d+2 



2(4vr)^ 



r 



d+l 



(if £(" 2 + 



N=l 



d+l 

7T 2 



+ j d+2 

r(-|) ^ 



(103) 



In order to compute the sum we use the Epstein-Hurwitz function which is defined as 

iv -r u, ; - h " 

AT=1 



( EH (s,a 2 ) = J2(N 2 + a 2 Y 



+ 



2 
2vr s 



2r( s ) 



[a' 2 



V)-f- 



T^n s "k s _i (^vmv 7 ^) , (104) 



n=l 



where is the modified Bessel function of second type. In our case s = — £ rr and a = -A 
It turns out that the second term in (11041) cancels with the second term in (11031) . Often this 
cancelation is not performed explicitly but by an equivalent argument the second term in 
the Epstein-Hurwitz function is discarded 41|. This is why people claim that in the zeta 
function regularization method it is not necessary to subtract any quantity and that a finite 
result comes directly considering only the first integral in ( )96l) . The final expression for the 
energy is 



d+2 



2(4^)* 



-ir 

2 



2 fl 2 x 1 



f- / +1 + -r V E -^K-mWn) • (105) 



The first term is a constant energy and therefore we discard it because does not contribute 
to the Casimir force. Deriving the energy with respect to the separation I between the 
hyperplanes, we obtain the (o?+2)-dimensional Casimir force 



d£, 



d+2 



d+2 



dl 



2 d+2 d 



/ A \ d + 2 



dl 



1 °° 1 

■j Yl -d±j K -^( 2nl f j ') 



h — i n 2 



(106) 



which can be evaluated explicitly. Using the properties K_ v (z) = K v {z) y zd z K v {z) 
—zK u _i(z) — uK u (z) we obtain finally 



f{d+2)(p) = -2 



47T 



1 + 2 



11 2 1 

J2 j+2 Kd±2 (2nln) --fV — Kd±4 (2n//i) 

1 1 ' <2 



d+2 

/ 2 — n 2 



n=l 



«2 



(107) 
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In particular we are interested in the 4D Casimir force, which is obtained setting d = 2 in 
the formula above 

,2 



fl- 



1 i o ^ i 

71=1 



n=l 



(108) 



Sometimes this expression is presented in a slightly different way, which can be obtained by 
using the identity K v (z) = K u _ 2 (z) + ~ K v -\{z) 



8tt 2 



oo _ _ OO -. 

n=l 



/ 2 ^— ' n : 

n=l 



(109) 



B. Green's function approach 



Our starting point is the integral ffTBj) . but in the spirit of generality, in analogy with the 
zeta function method we allow to have d transverse dimensions, namely 

7T 2 fc 2 

(110) 



1 f d d k 

}{d+2)\p) ~ 2 J (2^3 



9.-7T I / 



W^p 2 +/U 2 + 



7T 2 fc 2 
i 2 



This integral can be performed straightforward in polar coordinates, in terms of the volume 
of the unitary (^-dimensional sphere that we denote by vol(S d ) 



d d k f d£ vol(S d ) 



(2vr) d J 2n (2n) d + 1 J() 
In polar coordinates, Eq. (5) is rewritten as 
vol{S d 



p d dp, with, vol(S° 



in) 



/(<i+2) (aO 



fc=l 



2(27r) d+1 

vo /(^)r (*±I) 
4(27r) rf+1 

where we have used the result 



2 ^ ^2 ,00 



,2 1 .,2 1 7r 2 /c 2 ^ 



/ d PP d ^/p^T^ 

JO 



'1121 



E 

/c=i 



2Z.2\ — 



r (-**») 



<i+2 



„ *V + c) 



U/d + iX r( s -i-f) 



C 2 2 



(113) 



2 v 2 / r(s) 

with s = 1 and c = p 2 + Klf - in the first integral, and s = —1/2 and c = p 2 in the second. 
The next step in the computation is to notice that the two terms together in (11121) can be 
rewritten in terms of the derivative of a Epstein- Hurwitz function (11041) . In order to show 
this consider the infinite series in the first term which we shall denote as S 



1^7r 2 A; 2 /o n 2 k 2 \^ 1 /7r\^+id^/., l 2 p 2 ^^ 



I ^ I 2 

k=l 



P 2 + 



P 



d+l\l 



k=l 



;ii4) 
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but the series written in this way is precisely the Epstein-Hurwitz zeta function fjl04[) . In 
terms of it 



_ 1 {Tr\ d+1 d 

~J+l \l) Ji Ceh 

and computing explicitly the derivative we obtain 



d + 1 ( Ijj, 

71 



2 ' 



(115) 



S 



1 



d + ir(^±i) 



d 



20F ^ M ^ 2 ' , 



■ (H6) 



Using the identity -(d + l)T = 2T (-^f 1 ), and since T (-±) = -2^n, we obtain 

finally 



2ST 



+ tht dl 



oo 

d X \ d+2 



I 2 ^2 n 2 K_d+2 (2nlfM) 



;ii7) 



n=l 



But these are precisely the terms inside the brackets in equation (|112p . so we get the result 



l(S d )T(^)^ ^d 
(27i) d + 2 ^ 2 Jl 



oo 

d \_ — \ d+2 



d \ — * d+2 

I 2 2_^n 2 K d+2 (2n/ju) 



(118) 



n=l 



We can evaluate explicitly the derivative. Using the properties K- V (z) = K v (z) y 
zd z K v {z) = —zK v ^\(z) — vK v {z) and substituting the value of vol(S d ) we obtain finally 



/(d+2) (^) 



d+2 

/U \ — 

47T 



1 ^ 1 



Y] -4+2Kd+2(2nlfi) -|V -^jKd±4(2nlii) 
I 2 n 2 /2 n=1 n 2 



(119) 



This expression of the force coincides exactly with the one obtained above using the zeta 
function approach. 
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